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Ripplons - gravity- capillary waves on the free surface of a liquid or at the 
interfaces between two superfiuids - are the most favourable excitations for 
simulation of the general-relativistic effects related to horizons and ergore- 
gions. The white-hole horizon for the "relativistic" ripplons at the surface of 
the shallow liquid is easily simulated using the kitchen-bath hydraulic jump. 
The same white-hole horizon is observed in quantum liquid - superfiuid ^He. 
The ergoregion for the "non-relativistic" ripplons is generated in the exper- 
iments with two sliding ^He superfiuids. The common property experienced 
by all these ripplons is the Miles instability inside the ergoregion or horizon. 
Because of the universality of the Miles instability, one may expect that it 
could take place inside the horizon of the astrophysical black holes, if there is 
a preferred reference frame which comes from the trans- Planckian physics. 
If this is the case, the black hole would evapotate much faster than due to 
the Hawking radiation. Hawking radiation from the artificial black hole in 
terms of the quantum tunneling of phonons and ripplons is also discussed. 

PACS numbers: 04.70.-s, 67.40. Hf 



1. INTRODUCTION 

Black hole is the region of spacetime in which gravity is so strong that 
not even light can escape from this region. There are many astrophysical 
black hole candidates, but the unambiguous observational evidence for the 
existence of black holes has not yet been established (see, e.g. review papers 

The defining property of the black hole - the event horizon - is still 
missing in the experiment, though there are indications that some of the 
candidates might have event horizons. At the moment we are not able to 
study experimentally the exotic phenomena related to the event horizon. 
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for example, how the presence of the horizon modifies the properties of the 
vacuum of relativistic quantum fields. That is why the modeling of the event 
horizon in condensed matter systems may provide us with some ideas useful 
for study of the astronomical black holes, their formation and decay. The 
close analogy between the relativistic quantum vacuum and the ground state 
of quantum condensed matter appears to be instrumental. In this respect the 
quantum liquids, such as superfluids, are the best candidates for simulations: 
the horizon can be relatively easy constructed in the moving liquid, and the 
properties of the quantum vacuum in the presence of the horizon can be 
studied. 

Two-fluid hydrodynamics of superfluid liquids describes the coherent 
motion of the superfluid 'quantum vacuum' and dynamics of excitations 
above the ground state - quasiparticles - which form the normal component 
of the liquid. At low energy some quasiparticles mimic relativistic particles 
propagating in the effective curved space-time produced by the distortions 
of the superfluid vacuum and by its superfluid flow. This is an example of 
the quantum gravity provided by quantum liquids. Starting from the trans- 
Planckian physics - quantum many body system of interacting atoms - one 
obtains in the low energy corner the emergent quantum fields of relativistic 
quasiparticles living in the background of the emergent classical gravity field. 

The effective gravity in quantum liquids is dynamical, but as distinct 
from the Einstein theory of gravity, the metric field typically obeys the hy- 
drodynamic equations rather than Einstein equations. There is an exclusion 
from this rule: in superfluid '^He-A, the effective quantum electrodynam- 
ics of "relativistic" fermions and "photons" emerges at low energy. These 
quantum fields propagate in effective space-time whose metric field obeys 
equations remotely resembling Einstein equations^. All this occurs because 
superfluid ^He-A belongs to the special universality class of Fermi systems. 

However, even if the effective gravity in liquids is distinct from the 
Einstein gravity, it can be useful for simulations of many phenomena related 
to Einstein gravity, including the physics of black holes. For example, the 
effect of Hawking radiation is solely determined by the behavior of quantum 
fields in the vicinity of the event horizon, and it does not matter for which 
excitations (photons or sound waves) the horizon is constructed and from 
which equations (Einstein or hydrodynamic) it is obtained. 

The acoustic analog of a black hole - the horizon for sound waves - 
has been suggested by Unruh in 1981 The event horizon for sound waves 
emerges if the flow velocity of the liquid exceeds the local speed of sound 
c in some region of liquid: sound waves (or quasiparticles in superfluids - 
phonons) cannot escape from this region. The propagation of sound in the 
presence of the acoustic horizon is described by the so called acoustic metric, 
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which is similar to the metric describing the space-time in general relativity 
in the presence of the black hole. This is discussed in Sec. ID together with 
the analog of the Hawking radiation which is described in terms of quantum 
tunnelling between the classical trajectories'^^. Since 1981 there appeared 
many suggestions to simulate the black hole and white hole horizons for 
various excitations in various laboratory systems (see review paper ^ and 
references therein). 

In Sec. ED we discuss the most perspective analog - the effective 2+1 
space-time emerging for the surface waves (ripplons) in the shallow wa- 
ter limit. The effective horizon for ripplons has been first suggested by 
Schutzhold and Unruhl^. This analog demonstrates the mechanism of the 
decay of the black hole which is alternative to the Hawking radiation - the 
instability of the quantum vacuum behind the horizon (Sec. I3.4.() . This 
mechanism is applicable both to the relativistic systems with the horizon 
and non-relativistic systems with the ergoregion, and it has been experimen- 
tally observed for the non-relativistic ripplons propagating at the interface 
between two superfluids 1 9 1 ^9 1 ^ 3. Similar instability of the vacuum inside the 
astronomical black hole is possible. 

In Sec. Elwe discuss experiments with the hydraulic jump in superfluid 
^He^. This circular hydraulic jump simulates the 2+1 dimensional white 
hole for the relativistic ripplons and the instability inside the horizon 

In Sec. ED some perspectives are discussed. 



2. GRAVITY FOR PHONONS 

2.1. Effective Metric 

In the frame moving with velocity Vg of the superfluid vacuum - the 
comoving frame - the spectrum of sound waves is 'relativistic': u;™™(k) = 
±cfc, where c is the speed of sound. In the laboratory frame, the spectrum 
is Doppler shifted: 

a;(k,r) = ic/c + k- Vs(r) , (1) 

or 

(w - k • Vs)^ = c^fc^ . (2) 
This can be written in the general Lorentzian form 

g'^'^k^k, = , (3) 

where ko = —uj, and the contravariant components of the metric are 

^00 ^ _i ^ ^o,.(r) ^ , g'^ir) = cH'^ - vHMiv) . (4) 



G.E. Volovik 



The inverse matrix g^i, determines the trajectories of phonons in the moving 
hquid. This covariant metric describes the effective (3+l)-dimensional space- 
time in which phonons move along the geodesic hnes ds'^ = 0: 

ds"^ = g^udx^'dx'' = -dt^ + c"^(dr - Ws{r)dtf . (5) 



2.2. Black and White Holes 

The most instructive for us is the spherically symmetric flow when the 
superfluid velocity is radial: Vs(r) = rv{r): 

ds^ = -(1 - v'^{r)/c^)dt^ - 2c~^v{r)drdt + c-^idvf . (6) 

The coordinate transformation 

r v(r) 

t = t+ dr ^ ^ ^ . (7) 
leads to the more familiar metric 

ds' = -[l-^\d? + -^-—^ + -,r'd^' . (8) 

In the particular case when the velocity field has the form v'^{r) = 
2GM/r, the equation © reproduces the Schwarzschild line element in Ein- 
stein gravity, where M is the mass of the gravitating body, and G is the 
Newton constant. The corresponding equation Q with v'^{r) = 2GM /r is 
known as the Painleve-Gullstrand line element which describes the gravity 
field of the same body but in different coordinate system ~ the coordinate 
system adapted to the free falling particle 

In moving superfluids, the Painleve-Gullstrand type metric with 
nonzero off-diagonal element qqi is naturally generated by fiow since g^i oc v^i ■ 
The notion of the black and white event horizons has also a very simple 
meaning in liquids. For example, let us assume that we are able to construct 
the spherically symmetric flow with v{r) < such that within some radius 
r < rfi one has v^(r) > c^. Then phonons cannot escape from this region 
since they are dragged by the superfluid vacuum to the center with the speed 
faster than their velocity c. Thus the surface r = at which v{r) = —c 
would serve as the horizon of the black hole. 

On the contrary, if the liquid flows outward, i.e. for v{r) > 0, all the 
phonons are necessarily dragged away from the region r < r^. In this case, 
the surface r = r/j at which v{r) = c would serve as the horizon of the white 
hole. 
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Though in the Einstein gravity the Schwarzschild metric and the Painleve- 
Gullstrand metric describe the same gravitational field, the latter is more 
preferable especially when the physics of the event horizon is concerned. 
The Schwarzschild metric contains the singularity at the horizon: the q^t- 
element of the metric is infinite at f^(r) = (? . This is the so-called coor- 
dinate singularity - the non-physical singularity which can be removed by 
coordinate transformation to the Painleve-Gullstrand metric. The latter is 
smooth across the horizon, which allows us to study the behavior of the 
quantum fields and quantum vacuum both outside and inside the horizon. 
This is the reason why the Painleve-Gullstrand metric which provides us 
with the simple "river" model of black hole popular now in physics of the 
event horizon (see Ref . ^ and references therein) . 



As distinct form the Einstein equations which admit the spherically 
symmetric solutions with the horizon, the equations of hydrodynamics do 
not support the spherically symmetric acoustic horizon. Acoustic horizon 
occurs in a different geometry: in the so-called Laval nozzle, where the 
horizon takes place at the narrowest cross section of the tube. Fortunately, 
the shape of the surface of the horizon is not important for the discussion of 
the most interesting quantum properties of the quantum fields in the horizon 
- the Hawking radiation. This effect can be illustrated using the simplest 
case of the one-dimensional flow with varying velocity v{x) and varying speed 
of sound c{x). The effective acoustic metric experienced by phonons has the 
form 



and the horizon occurs at the point x = xq where v{xq) = c(xo), or v{xq) = 
— c(xo). The effective gravitational field at the horizon is given by 



Let us choose the flow with v{x) > 0. To simulate the black hole 
horizon, we must take v{x) < c{x) at x < xq (exterior region) and v{x) > 
c{x) aX X > xq (region behind the horizon). In the local reference frame 
comoving with liquid, the energy spectrum of phonons is everywhere positive 
TiLo'^°'^{k) = Tick > 0, and thus at T = phonons are not excited. This 
vacuum state is, however, not in equilibrium since in the local comoving 
frame the whole velocity field is time dependent and thus the quasiparticle 
energy and the vacuum as a whole is not well defined. The energy is well 



2.3. Hawking Radiation as Tunneling 






(10) 
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defined in the laboratory frame where the velocity field is time independent. 
But in this frame some modes acquire negative energy behind the horizon, 
^a;(k) = hck + /ik • Vg < 0. All this indicates that the initial comoving 
quantum vacuum is unstable in the presence of the black hole. Hawking 
radiation provides the mechanism for the decay of the black hole unless 
some other, more violent, process of vacuum instability intervenes. 

In the semiclassical description, the Hawking radiation can be consid- 
ered as the quantum tunneling between classical trajectories'^'^ (the same 
analysis applied for the Unruh effect in superfluids see in Ref. CEI). We 
consider the positive frequency modes, to = uj{k) > 0, as viewed in the lab- 
oratory frame. Classical trajectories k{x) of phonons may be found from 
Eq.(^: u) = ±c{x)\k{x)\ + v{x)k{x). There are two relevant trajectories in 
Fig. HI 

m = . / . . • (11) 

v[x) — c{x) 

The trajectory with k{x) < starts at the horizon and propagates into the 
exterior region. In the comoving frame the phonon frequency is positive, 
^com _ c(x)\k(x)\. Formally the trajectory starts with infinite frequency, 
a;™™(xo) = +00, but practically there is the cut-off which in quantum liquids 
is given by the Debye frequency, and in the quantum vacuum - by the Planck 
energy scale. The trajectory with k{x) > also starts near horizon but 
propagates into the region behind the horizon. In the comoving frame the 
phonon frequency is negative, uj^°^{k) = —c{x)\k{x)\, and approaches — cxd 
at the horizon. 

Quantum mechanics allows the phonon to tunnel between the two tra- 
jectories. The tunneling exponent is determined by the conventional semi- 
classical action 

S = lu.[nkix)dx = lm[^^^dx, (12) 
J J v{x) — c[x) 

The momentum k{x) on the classical trajectory has a pole at the horizon 

„ _Jp^ , ,13) 

gx{x - xoj 

Shifting the contour of integration into the upper half-plane of the complex 
variable x, one obtains for the tunneling action 

^^TThucixol ^^^^ 

9x 

The tunneling between the trajectories describes the process of creation of 
the pair of phonons, one of which is radiated into the exterior region. The 
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v(x) < c(x) 



v(x) > c(x) 



X 



Fig. 1. Two relevant rajectories of massless 'relativistic' quasiparticles - 
phonons - travelling from the horizon of acoustic black hole. Arrows show 
the direction of motion. Quasiparticles have diverging energy in the vicinity 
of the horizon (infinite blue shift) which must be cut off by Planckian (De- 
bye) energy scale. The tunneling between the two trajectories describes the 
process of emission of phonons from the horizon which is analogous to the 
Hawking radiation from the black hole. 

probability of tunneling demonstrates that the radiation from the (acoustic) 
black hole is thermal with the Hawking temperature Th being proportional 
to the (effective) gravity at the horizon: 



This semiclassical approach is only valid when the tunneling action is large 
compared to Planck constant, i.e. when huj ^ Th- 

The Hawking radiation is completely determined by the (effective) met- 
ric in the vicinity of the horizon. It does not depend on the motion equations 
for the gravity field: the metric can obey Einstein, hydrodynamic or other 
equations which allow for the existence of an event horizon. The Ea. (|15() is 
applicable for any relativistic excitations radiated from the horizon, includ- 
ing 'relativistic' fermionic quasiparticles living in superfluid ^He-A EEI_ 

Note that in the above derivation of the Hawking temperature we ac- 
tually did not use the full quantum mechanics. Only the wave mechanics 
was exploited, which is also applicable to the classical waves propagating in 





(15) 
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the background classical metric. Ea. HlSjl can be rewritten without using the 
Planck constant h (see Ref. 

Pocexpi — —] , a^H = TT^ • (16) 

The above equation describes the spectrum of classical waves radiated from 
the horizon. In the full quantum mechanics the radiation is caused by quan- 
tum fluctuations in the vicinity of the horizon. But in the wave mechanics 
the radiation can be generated by some external noise introduced in the 
vicinity of the horizon or by thermodynamic fluctuations. However, these 
initial perturbations near the horizon are only responsible for the prefactor 
in Ea. H16() . which was missing in our derivation, but they do not influence 
the value of the Hawking frequency wh- That is why the Hawking frequency 
can be measured in classical experiments. 



2.4. Does Hawking Radiation Exist? 

Let us mention that in the literature there is a controversy concerning 
the Hawking radiation. There is even a point of view that a Schwarzschild 
black hole formed during a collapse process does not radiate (see e.g. Ref. 

The other authors (see e.g. Ref. suggest that the Hawking radiation 
may exist but the Hawking temperature is twice larger than that originally 
calculated by Hawking, i.e. twice that in Ea.((T3|). 

Prom our point of view the existence of Hawking radiation and the 
Hawking temperature depend on the vacuum state. In our case when the 
horizon is formed by the moving superfluid liquid, the initial vacuum is 
the comoving vacuum, i.e. the vacuum seen in the frame moving with the 
liquid. The comoving vacuum does not coincide with the equilibrium vacuum 
determined by the Killing vector energy (the vacuum which is seen in the 
laboratory frame). That is why the comoving vacuum is unstable, and the 
initial stage of the development of instability is the Hawking radiation (see 
discussion in Ref. These two vacua and the quantum tunneling are 

described in the Painleve-Gullstrand frame, which is the natural frame for 
the propagation of quasiparticles in the moving liquids. 

On the contrary, the Schwarzschild coordinates are irrelevant for this 
problem: they are obtained after singular coordinate transformation which 
is unphysical in superfluids since it removes the part of the physical space. 
This is the typical situation in condensed matter, that in general the singular 
coordinate transformations are forbidden. As a result the coordinate frames 
are distributed into classes: the transition between the classes is only pos- 
sible after singular ccordinate transformations. Por the artificial (acoustic) 
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black holes where the proper frame is the Painleve-Gullstrand frame and 
the proper initial vacuum is the comoving vacuum, the Hawking radiation 
is the real process. The original Hawking's derivation has been made in 
the frame which belongs to the same class as the Painleve-Gullstrand frame. 
That is why the radiation from the artificial black hole is characterized by 
the same temperature as was derived by Hawking. However, the Hawking 
radiation from the astrophysical black hole stiil remains the open question 
since the problem of the correct choice of the vacuum state after collapse is 
not solved. 

Note also that for quantum liquids, in some cases the Hawking radiation 
is suppressed because of the effect which is similar to the Coulomb block- 
ade in mesoscopic physics: the radiated quasiparticle disturbs the vacuum 
state so that the radiation becomes impossible. For example, the radiated 
quasiparticle changes the superfluid current and thus violates the mass con- 
servation law. However the radiation of two quasiparticles is possible, since 
it does not violate the conservation law. The simultaneous tunneling of two 
quasiparticles is called the co-tunneling. The corresponding tunneling expo- 
nent describing the co-tunneling process is twice the tunneling exponent of 
a single quasiparticle. As a result the Hawking temperature is twice smaller 
than in Ea. ((T3|l . i.e the Coulomb blockade effect is opposite to that suggested 
in Ref. m 



3. GRAVITY FOR RIPPLONS 

In all the condensed matter systems suggested for simulation of gravity, 
even if the horizon can be constructed in principle, the estimate for the 
Hawking temperature is very pessimistic for the observation of the Hawking 
radiation (though the classical analog of the Hawking effect discussed in 
the previous section can be observed). That is why it is instructive to look 
for the other possible experimental consequences of the event horizon. The 
effective gravity for ripplons - the capillary-gravity waves on the free surface 
of the liquid or at the interface between two liquids - provides us with the 
new mechanism of the decay of the black hole. Let us start with ripplons 
propagating on the surface of the liquid. 

3.1. Effective Metric on the Surface of Liquid 

The general dispersion relation <-ij(k) for ripplons on the surface of a 
liquid is 

M{k){uj -k■^^f = F + k'^a . (17) 
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Here a is the surface tension; F = pg is the gravity force where p is mass 
density of the hquid; and v is the velocity of the hquid. The quantity M{k) 
is the /c-dependent mass of the hquid which is forced into motion by the 
osciUating surface: 

= triirrk • ^''^ 

where h is the thicknesses of the layer of the liquid. 

The spectrum H17|) becomes "relativistic" in the shallow water limit 
kh<^l, k<^ ko: 

{u;-h■^rf=c^k' + c^k'(^■^^-h'^ , c' = gh , kl = pg/a. (19) 

If the k'^ corrections are ignored, the spectrum of ripplons is described by 
the effective 2+1 metric l^l 

g^'^'k^k^ = , k^, = {-uj, k^, ky) , (20) 

with the following elements 

^00 = -1 , = -yi , gij = c^S'^ - v'v^ . (21) 

The interval describing the effective 2+1 space-time in which ripplons prop- 
agate along geodesies and the corresponding covariant components of the 
effective metric are 

ds'^ = g^vdx^dx" , g^Q = -I + — , g^^ = -— ^ g^j = —5ij . (22) 



3.2. Effective Metric on the Interface Between Two Superfluids 

The spectrum of ripplons propagating along the interface between two 
superfluids: 

Mi{k){uj - k • vi)2 + M2{k){uj - k • va)^ = F + k^a . (23) 

Here, as before, a is the surface tension of the interface and F is the force 
stabilizing the position of the interface. For the interface between ^He-A 
and ^He-B the gravity force is negligible, since these liquids have almost 
equal densities: \pi — /O2I ~ 10~^(/3i + p2)- However, these two superfluids 
have essentially different magnetic properties, and the corresponding force 
F which stabilizes the position of the interface is provided by the gradient 
of external magnetic field. The effective masses Mi (A;) and M2{k) are : 
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where hi and /i2 are the thicknesses of the layers of two superfluids; psi 
and ps2 are superfluid densities of the liquids. Further we assume that the 
temperature is low enough so that the normal fraction of each of the two 
superfluid liquids is small, then psi ~ Pi and ps2 ~ P2- 

Experiments with the AB interface 1^ where conducted in the 'deep wa- 
ter' limit khi ^> 1 and kh2 ^ 1. In the opposite limit of a thin slab, where 
khi <C 1 and kh2 <^ 1, one obtains 

ai(u; - k • vi)2 + ^2(0; - k • va)^ = c^A;^ , (25) 

where 

h2Pi hip2 2 Fhih2 

"1 = 1 ' "2 = T -j^ , c = — . 26 

h2Pi + hip2 h2Pi + hip2 h2Pi + hip2 

This can be rewritten in the Lorentzian form (|2()j) with the following effective 
contravariant metric g^^: 

gOO = _i , gOi = _aiv{ - a2v\ , g^^ = c^6'^ - aiv\v{ - a2v\v{ . (27) 



3.3. Interaction with Environment 

The spectra (|T7j) and are valid for the perfect fluid, where dissi- 
pation due to friction and viscosity is neglected. They must be modified 
when the dissipation is added. For the ripplons propagating at the interface 
between two superfluids t he di ssipation leads to a simple extra term on the 
right-hand side of Ea.(l23l)M 

Mi{k)[uj - k • vi)2 + M2{k){L0 - k • ^2? = F + k^a- iVoj . (28) 

For the ripplons at the interface between ^He-A and ^He-B the friction pa- 
rameter r > depends on temperature and is proportional to at low T^^. 
The important property of the added dissipative term is that it introduces 
the reference frame of the environment. The w-dependence of the dissipative 
term in Eq. (|28|) . which has no Doppler shift, implies that this spectrum is 
written in the frame of the container. 

Since the main effect of this term is the introduction of the distinguished 
reference frame, the explicit form of this term is not important. That is why 
we expect that similar dissipative interaction with the environment exists 
for ripplons propagating on the free surface of the superfluid and even on 
the surface of normal viscous liquid: 

M{k){oj -^k■^^f = F + k'^a -iVu . (29) 

The parameter P can be considered as phenomenological parameter, which 
depends on cj, k and the Reynolds number of the flowing liquid. 
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3.4. Instability in the Ergoregion and Landau Criterion 

In superfluids the stationary (time- independent) flow of superfluid com- 
ponent is frictionless until the critical velocity of flow is reached. This is the 
essence of the phenomenon of superfluidity. However, the notion of criti- 
cal velocity implies that there is a preferred reference frame with respect 
to which the velocity is counted. This is the reference frame of the envi- 
ronment (the frame of the container). In equations, this reference frame is 
provided by the F-term, which describes the dissipative interaction with the 
environment. This term also gives rise to the attenuation of ripplons. 

When this term is taken into account, from the spectrum uj{k) in Eq. 
()29() it follows that the instability to the formation of the surface waves 
occurs when the velocity v of the flow with respect to the wall exceeds 
the critical velocity vl, at which the imaginary part Im u{kc) of the energy 
spectrum of the critical ripplon with momentum kc crosses zero and becomes 
positive aX V > Vl. At that moment the attenuation of ripplons transforms 
to amplification; and critical ripplons start to grow exponentially 1^^. 

The critical velocity vl and the momentum of the critical ripplon kc do 
not depend on the friction parameter T. Moreover, at this velocity the real 
part Re u;(k) = a;™™(/c) -|- k • Vg also crosses zero: it becomes negative at 
V > Vl (see Fig. |2lfor ripplons at the AB interface, which we will discuss in 
Sec l3.5.() . This means that the threshold velocity coincides with the Landau 
criterion for the ripplon nucleation: 

VL = miufe— , ^com(^) ^ ^^pg ^ ak^)/M{k) . (30) 

The Landau critical velocity is different in the "relativistic" shallow-water 
and "non-relativistic" deep-water regimes: 

Vl = c , kc = , if hko < , (31) 
vl = c\j2/hkQ , kc = ko , if hko » 1 . (32) 

In both regimes the frequency of the critical ripplon is Lo{kc) = 0, i.e. the 
critical ripplon must be stationary in the wall frame. 

The region in the liquid, where the flow velocity v exceeds vl, represents 
the ergoregion, since in the container frame the energy of the critical ripplon 
is negative in this region, a;™™(A;) -|-k- v < 0. For the "relativistic" ripplons, 
the ergoregion is the region where v exceeds c. This can be expressed in 
terms of the effective metric in Eq. (|22|) : in the ergoregion the metric element 
goQ changes sign and becomes positive. If the flow is perpendicular to the 
ergosurface (the boundary of the ergoregion), then the ergosurface serves 
as the event horizon for ripplons. It is the black hole horizon, if the liquid 
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moves into the ergoregion, and it is the white hole horizon if the hquid 
moves from the ergoregion. However, the exponential growth of ripplons in 
the supercritical region is universal and does not depend on direction of the 
flow in the vicinity of the ergosurface. The flow can represent the white-hole 
or black-hole horizons; the relativistic ergoregion without the horizon; and 
even the ergogerion for the non-relativisic quasiparticles. 

The instability develops due to the interaction with the fixed reference 
frame and occurs in the region where the energy of the critical fluctuation is 
negative in this frame. Such kind of instability is also called the Miles insta- 
bility. I^In principle, Miles instability may take place behind the horizon of 
the atsronomical black holes if there exists the fundamental reference frame 
related for example with Planck physics EEH. It may lead to the decay of the 
black hole much faster than the decay due to Hawking radiation. 



3.5. Ergoregion Instability vs Kelvin-Helmholtz Instability 

In the case of the interface between ^He-A and ^He-B, the critical veloc- 
ity of instability towards the growth of critical ripplon has been measured 
in the nonrelativistic deep-water regime ^ and has been found in a good 
agreement with the theoretical estmate of the Landau velocity (modified for 
the case of two liquids I^El) wihout any fitting parameter (see Ref. where 
also the turbulence triggered in bulk '^He-B by the interface instabilty is 
discussed). The important additional feature, which takes place in case of 
the interface between two sliding superfluids, is the existence of two different 
critical velocities. 

One of them is the Landau critical velocity, above which the ergoregion 
appears for ripplons. In the typical situation the A-phase is stationary with 
respect to the walls of container, while the B-phase is moving with respect 
to the container with the velocity VsB- For the deep-water regime and in the 
limit of low temperature T <C Tc, this critical velocity is given by: 

vl = -^Vf^ . (33) 

The second critical velocity occurs if one completely neglects the interaction 
with the enviroment putting the friction parameter F = 0. In this case 
the preferred reference frame disappears from the equations, and thus the 
physics is determined by the relative velocity vi — V2 of two superfluids. 
In this case it is well known that the Kelvin-Helmholtz instability of the 
interface occurs when the velocity difference reaches some threshold value. 
In our case, when the A-phase is stationary with respect in the container 
frame, while the B-phase is moving with respect to the container with the 
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slope 




PVsB =4V Fa 



Landau criterion 
(ergoregion instability) 



Kelvin-Helmlioltz criterion 
(physical singularity) 



Fig. 2. Two critical velocities for ripplons living at the interface between two 
sliding superfluids - A and B phases of '^He. The lowest one is the Landau 
critical velocity, at which the energy of ripplons first becomes negative. The 
region with negative energy is called the ergoregion. If the friction parameter 
r describing interaction with the environment is non-zero, the imaginary 
part of the frequency becomes positive in the ergoregion, i.e. attenuation 
of ripplons transforms to the amplification in the ergoregion which indicates 
the instability with increment proportional to T. The second velocity marks 
the more crucial Kelvin-Helmholtz instability at which the imarginary parts 
of both signs emerge. In the 'relativistic' shallow-water regime this velocity 
corresponds to physical singularity in metric field (Fig. ISJ- The shown 
values of two critical velocities correspond T <C T^, deep-water limit, and 
the experimental condition for the flow of two superfluids: the A-phase is 
stationary with respect to the walls of container, and the B-phase is moving 
with respect to the container with the velocity VsB- 

velocity VsB, the relative velocity of two superfluids is |vi — V2I = VsB and 
the Kelvin-Helmholtz criterion becomes 



Here we took into account that the densities of two liquids are equal pi = 
P2 = p. The difference between two critical velocities in such an arrangement 



The behavior of the real and imaginary parts of the spectrum of the 
crtitical ripplon as functions of VsB is shown in Fig. [21 When F 7^ 0, at vkh 
one has the smooth crossover, which becomes the true threshold at F = 0. 
In the "relativistic" shallow-water regime, both thresholds have physical 
meaning, especially when the flow is radial (Fig. The circumference 
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Fig. 3. Artificial 2+1 black hole for ripplons at the AB interface in the 
"relativistic" shallow-water limit, which can occur for the radial flow of the 
B-phase near the interface. The black hole horizon occurs when the flow 
reaches the Landau critical velocity \v{r)\ = vl- At larger critical velocity 
\v{r)\ = vkh > Vl the Kelvin-Helmholtz instability occurs, which in terms 
of effective metric corresponds to the physical singularity within the black 
hole. 



v{r) = —Vl = — c/-y/ai represents the black-hole horizon. Deep inside the 
black hole, at v{r) = —vkh = — c/yaia2, the real singularity takes place: 
the g^r and gor components become infinite together with the determinant 
g of the metric g^y. 



4. HYDRAULIC JUMP AS WHITE HOLE 

4.1. Hydraulic Jump 

The analog of a white hole horizon for the 'relativistic' ripplons on the 
surface of a shallow liquid is achieved in the kitchen-bath experiments. This 
is the so-called hydraulic jump first discussed by Rayleigh in terms of the 
shock wavel^. The circular hydraulic jump occurs when the vertical jet of 
water falls on a flat horizontal surface (Fig. IIJ. The flow of the liquid at 
the surface exhibits a ring discontinuity at a certain distance r = R from 
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Fig. 4. Hydraulic jump. The standing waves at the otherwise smooth surface 
inside the white hole is the result of the egroregion instability. Courtesy Piotr 
Pieranski 

the jet (observation of the non-circular hydraulic jumps with sharp corners 
has been reported in Ref. At r = R there is an abrupt increase in the 
depth h of the liquid (typically by order of magnitude) and correspondingly 
a decrease in the radial velocity of the liquid. The velocity of the liquid in the 
interior region (r < R) exceeds the speed of 'light' for ripplons v > c = \/Tig, 
while outside the hydraulic jump (r > R) one has v < c = ^/hg (Fig. [SJ- 
Since the velocity flow is radial and outward, the interior region imitates the 
'white-hole' region. The interval of the 2+1 dimensional effective space-time 
in which the "relativistic" ripplons "live" is 

ds^ = -c^dt^ + {dr - v{r)dtf + r'^d<i)'^ . (35) 

The surface inside the white hole is smooth, since all the perturbations 
are removed by 'superluminal' flow of the liquid. In general relativity the 
metric is continuous across the horizon. In our case there is a real physical 
singularity at the white-hole horizon - the jump in the effective metric (|22|) . 
However, the discussed mechanism of the Miles instability in the ergoregion 
(or behind the horizon) discussed in Sec. IH.4.1 towards the growth of the 
critical ripplon does not depend on whether the horizon or ergosurface is 
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smooth or singular. According to this instabihty the critical ripplon must 
be stationary in the wall frame. The standing waves inside the white hole 
are clearly seen in Fig. HI The exponential growth of the critical ripplon 
is saturated by the non-linear effects, and then the whole pattern remains 
stationary (time- independent but not static). 

A similar 3D analog of the black or white hole with the physical sin- 
gularity at the horizon has been discussed by Vachaspati The role of 
the physical horizon is played by the phase boundary between two quantum 
vacua - two superfluid phases of the same liquid, such as ^He-A and ^He- 
B. Two vacua have different 'speed of light', say, different speed of sound, 
ci 7^ C2. The analog of black or white horizon occurs if the superfluid velocity 
of the flow through the phase boundary is subsonic in one of the superfluids 
but supersonic in the other one, ci < v < C2- 



4.2. Hydraulic Jump in Superfluids 

The analogy between the instability of the surface inside the hydraulic 
jump and the instability of the vacuum behind the horizon can be useful 
only if the liquid simulates the quantum vacuum. For that, the liquid must 
be quantum, and its flow should not exhibit any friction in the absence of 
a horizon. That is why the full analogy could occur if one uses either the 
flow of quantum liquid with high Reynolds number, or the superfluid liquid 
which has no viscosity. Quantum liquids such as superfluid or normal ^He 
and ^He are good candidates. 

The first observation of the circular hydraulic jump in superfluid liquid 
(superfluid ^He) was reported in Ref. The surface waves generated in 
the ergoregion (in the region inside the jump in Fig. [SJ were also observed. 
As in the normal liquid, this instability is saturated by the non-linear terms. 
The same instability observed at the interface between A-phase and B-phase 
has a different consequence the instability is not saturated and leads to 
the crucial rearrangement of the vacuum state: quantized vortices start to 
penetrate into the ^He-B side from ^He-A^Sl. They partially or fully screen 
the ^He-B flow and reduce its velocity back below the threshold for the 
ripplon formation. Thus the instability kills the ergoregion. 

Under the conditions of experiment ^ the hydraulic jump in superfluid 
^He is very similar to that in the normal liquid ^He. The position R of the 
hydraulic jump as a function of temperature does not experience disconti- 
nuity at the superfluid transition. This suggests that quantized vortices are 
present in the flow providing the mutual friction between the superfluid and 
normal components of the liquid. As a result even below the A-point, the 
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Fig. 5. Instability in the ergoregion. At the horizon both the energy and 
imaginary part of the ripplon spectrum cross zero. In the ergoregion where 
the energy is negative, the imaginary part is positive indicating instability 
towards the growth of the critical ripplon. 



liquid moves as a whole though with lower viscosity because of the reduced 
fraction of the normal component. 

To avoid the effect of the normal component it would be desirable to 
reduce the temperature or to conduct similar experiments in a shallow su- 
per fluid ^He. 

The advantage of superfluid ^He is that, as distinct from the superfluid 
^He, vortices are not easily formed there: the energy barrier for vortex nucle- 
ation in ^He-B is about 10^ times bigger than temperature In addition, 
in superfluid ^He the normal component of the liquid is very viscous com- 
pared to that in superfluid ^He. In the normal state the kinematic viscosity 
is v ^ 10~^ cm^/s in liquid ^He, and ~ 1 cm^/s in liquid ^He. That is 
why in many practical arrangements the normal component in superfluid 
^He remains at rest with respect to the reference frame of the wall and thus 
does not produce any dissipation if the flow of the superfluid component is 
sub-critical. 

One can also exploit thin films of a superfluid liquid, where the normal 
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component is fixed. The ripplons there represent the so-cahed third sound 
(recent discussion on the third sound propagating in superfluid ^He films 
can be found in Ref. In 1999 Seamus Davis suggested to use the third 
sound in superfluid ^He for simulation of the horizons 

In normal liquids it is the viscosity which determines the position R of 
the hydraulic jump (see^. The open question is what is the dissipation 
mechanism which determines the position R of the white-hole horizon in a 
superfluid flow with stationary or absent normal component when its viscos- 
ity is effectively switched off. Since there is no dissipation of the superfluid 
flow if its velocity is below vl, one may expect that the same mechanism, 
which is responsible for dissipation in the presence of the horizon, also deter- 
mines the position R of the horizon. If so, the measurement of R as function 
of parameters of the system will give the information on various mechanisms 
of decay of white hole. It is also unclear whether it is possible to approach 
the limit of a smooth horizon, without the shock wave of the hydraulic jump; 
and whether it is possible to construct the inward flow of the liquid which 
would serve as analog of the black hole horizon. 



5. CONCLUSION 

It appears that the ripplons - gravity-capillary waves on the surface of 
liquids or at the interfaces between superfluids - are the most favourable 
excitations for simulation of the effects related to horizons and ergoregions. 
The white-hole horizon for the "relativistic" ripplons at the surface of the 
shallow liquid is easily simulated using the kitchen-bath hydraulic jump. 
The same white-hole horizon is observed in quantum liquid - superfluid ^He 

The ergoregion for the "non-relativistic" ripplons is generated in the 
experiments with two sliding superfluids'^. 

The common property experienced by the ripplons in all these cases is 
the Miles instability inside the ergoregion or horizon. In some cases these 
instability is saturated leading to the standing waves inside the horizon or 
ergoregion, while in the other cases, such as for ripplons at the interface 
between ^He-A and ^He-B, the instability leads to the complete elimination 
of the ergoregion. Because of the universality of the Miles instability, one 
may expect that it could take place inside the horizon of the astrophysical 
black holes, if there is a preferred reference frame which comes from the 
trans-Planckian physics (see Sec. 32.4 in Ref.^. If this is the case, the black 
hole would evapotate much faster than due to the Hawking radiation. 

I hope, the future experiments in quantum liquids will explore the quan- 
tum limit, where the quantum effects related to horizon and ergoregion be- 
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come more pronounced. Experiments with hydraulic jump in superfluid He 
must be extended to the low-temperature region T <C Tc, while the ex- 
periments with ripplons at the interface between ^He superfluids must be 
extended to the 'shallow- water' regime, where ripplons become 'relativistic'. 
And, of course, ripplons at the interface between different superfluid states 
of ultra-cold Bose and Fermi gases shall be used. 

If the vacuum instability inside the horizon/ergoregion is saturated as 
in experiments with a free surface the other mechanisms will intervene 
such as the black-hole laser and even the quantum mechanical Hawking 
radiation of ripplons. The latter should be enhanced at the sharp discontin- 
uous horizon of the hydraulic jump and maybe near the sharp corners of the 
non-circular (polygonal) hydraulic jump observed in Ref. 
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